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Abstract 

We study asymptotic distribution of eigen- values to of a quadratic operator poly- 
nomial of the following form {uj 2 — L(uj))(pu = 0, where L(uj) is a second order differ- 
ential positive elliptic operator with quadratic dependence on the spectral parameter 
uj. We derive asymptotics of the spectral density in this problem and show how to 
compute coefficients of its asymptotic expansion from coefficients of the asymptotic 
expansion of the trace of the heat kernel of L(u>). The leading term in the spectral 
asymptotics is the same as for a Laplacian in a cavity. The results have a number of 
physical applications. We illustrate them by examples of field equations in external 
stationary gravitational and gauge backgrounds. 
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1 Introduction 



This work is motivated by studying quantum effects in stationary gravitational or gauge 
classical background fields 0,0. If quantum fields <p are free or one is restricted by 
one-loop approximation the fields obey second order differential equations. In the cor- 
responding coordinates (x°,x % ) these equations can be brought to the following general 
form: 

~d 2 + L(ido)] = 0. (1.1) 

L(ido) is the differential operator which includes both time derivative do = d/dx° and 
space derivatives, 

L(id ) = L 2 + (ido)Li + (id ) 2 L , (1.2) 

Lk are fcth order differential operators with the space derivatives only. This structure of 
L(ido) can be inferred by analyzing dimensionalities. In what follows we assume that Lk 
act over a (^-dimensional compact space Aid- The concrete form of ( |1 . 1| ) will be specified 
in next section. 

In time-independent backgrounds, field excitations w (x°,af) = e~ wx ° 4>u(x l ) with en- 
ergy u) are solutions to 

u 2 - L(u)] <f) w = 0, (1.3) 

L(u) = L 2 + ujL 1 + uj 2 L q . (1.4) 

Equation (|1~3| ) is an eigen-value problem on the spectrum of to. If L x and L in ( |I~2| ) 
are absent one has a standard problem of finding eigen-values A = u 2 of a second order 
differential operator L 2 . In general, however, operators L\ and Lq in ( |1.3[ ) are non-trivial 
and all three operators Lk do not commute. Thus, ( |1.3|) depends polynomially on the 
spectral parameter u. 

In this paper we call ( |1.3|) the non-linear spectral problem (NLSP) (this should not 
be confused with problems which depend non-linearly on eigen- functions) . Equations 
like (|I~3|), ( |1.4| ) belong to the spectral theory of polynomial operator pencils [|J, a field 



of mathematics where important pioneering results were established by Keldysh Q fifty 
years ago. In addition to already mentioned applications, quadratic and more general 
operator polynomials appear in other physical problems, for instance, oscillations of a 
viscous fluid, Schrodinger equation with energy-dependent potential and etc. 

The aim of this work is to analyze asymptotic distribution of the spectrum of ( |1.3|) at 
large uj. This can be done by studying asymptotics of the spectral density. If the spectrum 
is real the information about its properties can be also derived from the function 

K(t) = 1 -J2e- t -\ f>0, (1.5) 



where ui are eigen-values of (|1.3j ). The behavior of K(t) at small t is connected with the 
distribution of large uj. In next sections we show that at small t, if L{oj) is a positive 
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elliptic operator, 

1 oo 

^(^T^Ek + fi/^], (1.6) 

This expansion has the same form as asymptotics of the trace of the heat kernel of a second 
order differential operator. Moreover, the coefficient ao is proportional to the volume of 
Aid while other a n and b n are local functionals of the external background fields and 
can be computed by using heat-kernel coefficients of the operator L{uj) in (|1.3|). Thus, 
the behavior of the spectrum of the NLSP at large uj 2 is very similar to behavior of the 
spectrum of a Laplacian in a cavity. 



There are several important applications of ( |1.6| ) in quantum field theory. For instance, 
the spectrum of uj defines the vacuum energy E and, in case of a finite temperature 
the free-energy of system F(f3). For Bose statistics 

E = W" , F(0)=/T 1 5>(l- e -*'). (1.7) 

Asymptotics (|1.6|) determine the form of divergence of E at large uj. One can also use 
( |1.6|) to derive the behaviour of F(j3) at large temperatures. For further discussion see 



There are two main approaches to the operator polynomials discussed in the literature 
(see H): the method of linearization and the method of factorization. The method of 
linearization reduces the spectral problem like (|1.3|) to a standard linear problem but in 
an extended Hilbert space. This method is analogous to the reduction of an nth order 
differential equation to a system of n first order equations. The factorization method 
reduces an operator polynomial to a product of pseudo-differential operators each of which 
depends linearly on the spectral parameter. 

The method used in this paper is different. Our idea is to find the spectrum of (|1.3| ) 
from the spectrum of L(uj), by considering at first uj as an independent parameter. This 
way seems to be more simple (at least for discussed class of polynomials) and it enables 
us to use the results of the spectral theory of elliptic operators. Another new element 
which is absent in the Keldysh approach is in using the inner product <, > analogous 
to the Klein-Gordon product known in relativistic field theory. This step is motivated 
by physical applications and is important for our analysis. The eigen-functions (j)^ with 
different uj are orthogonal with respect to <, > and the set of eigen-functions can be 
orthonormalized. 

Asymptotics ( |1.6|) were first reported in our earlier publication || for particular prob- 
lem related to fields in stationary but not static backgrounds. The aim of this work is to 



make arguments leading to ( |1.6|) more general. 

The paper is organized as follows. In next section we give more detailed definition of 
the class of NLSP which will be considered here. Conditions which restrict this class are 
dictated by physical requirements to equation (|1 . 1| ) . We then derive a relation between 
the spectral density of eigen- values of (|1.3|) and the spectral density of operators L(uj). 
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This relation is used in section 3 in studying spectral asymptotics of the NLSP and for 
calculation of coefficients a n and b n in ( |1.6|) . One of the consequences of these results 
is that for certain n, depending on the dimension d, coefficients a n , b n are related to 
heat kernel coefficients of the operator 8$ + L(id ), see section 4. Section 5 is devoted to 
examples of NLSP where the spectrum has a simple form which allows one a direct check 
of the obtained formulas. One example is a charged field in a constant gauge potential, 
the other is a field in a rotating Einstein universe R 1 x S 3 . Finally, section 6 is devoted 
to spectral problems related to the Dirac equation. Although our results are applicable 
to these problems additional consideration is needed because the corresponding operator 
L(u>) is not self-adjoint. 



2 Formulation of the problem 

We consider free fields in a stationary space-time of a dimension D — d + 1. Equation 
( |1 . 1|) for is reduced to problem (|1.3| ), (|1.4| ) where L(u) will be taken in the following 
form 

L(u) = -(V fc + iA k + iua k ){V k + %A k + iuja k ) +uB + V. (2.1) 

It will be assumed that L(uj) acts on sections to some vector bundle over Aid where Aid is 
a compact <i-dimensional Riemannian manifold with metric h^. Index k in (|2.1|) is raised 
and lowered by using the metric tensor hik, are the corresponding connections on Aid 
and Ak, ctfc, B and V are some matrix- valued fields. 

In what follows we assume that the spectrum of (|1.3| ) is real and consider w as a 



real parameter. We also make two additional assumptions: i) L{u) is a positive elliptic 
operator which is self-adjoint on a Hilbert space L 2 (.Md) on Aid and L 2 = L(0) has 
strictly positive spectrum; ii) for solutions of (|1.1| ) there can be defined a time-independent 
product ( , ) such that all positive frequency eigen-functions a; (x o ,£*) = e~ luJX ° (j) UJ {x % ) 
(uj > 0) have positive norm ((<f> w , > 0) while the norm of negative-frequency functions 
(uj < 0) is negative. According to (i) there are no zero modes (u = 0) in the spectrum. 

Positivity of £(u;) implies that Aid is a Euclidean manifold. The second assumption is 
dictated by physical considerations. To quantize the theory the field is divided onto two 
parts. The part which has a positive norm is connected with particles while the negative 
norm part corresponds to anti-particles. Our assumption guarantees that energies uj of 
all physical particles are non-negative. 

The Hilbert space L 2 (A^d) is equipped with the standard inner product (0, tp) = 
f 4> + ip\/hd d x. The product ( , ) is defined as follows 

(0, ij>) = z(0, - i{<j>, fj)) - (0, LiV) - L Q i>) + z(0, LqV), (2.2) 



where = <9o0 and Lk are Hermitian operators defined for by Q1.4Q . One can check 
that for any two solutions and ip of which belong to L 2 (A^d) (|2.2j ) does not depend 



4 



on x°, do((j), ip) = 0. In a covariant theory on a globally hyperbolic space-time the latter 
property is extended to independence on the choice of a space-like hyper-surface. For 
scalar fields (|2.2|) coincides with the Klein-Gordon product. 
Consider the eigen-value problem for operator L(uj) 

L(cu)^ = A(w)^\ (2.3) 

where uj is a real parameter. If the spectrum A (a;) is known the spectrum of NLSP (|1.3| ) 
is determined by roots of equation 

X(",A) = 0, (2.4) 

X( U ,A) = W 2 -A( W ). (2.5) 
This also yields the eigen functions of the NLSP, W = 0^4 'a ; where uj is a root of 



and C w is a normalization coefficient. (Recall that ( |2.4j ) is assumed to have only real 
roots.) It is not difficult to see that for eigen-functions of (|1.3|) the following property 
takes place 

(<f>u>,i>o) = <L<7 x'(^)(0o;,^), (2.6) 

where 5^ = if uj ^ a and <5 W(7 = 1 if uj = o. x\uj) denotes the derivative 9 w x(o;,A) 
computed at a; which are roots of ( |2.4| ). It follows from ( |2.6| ) that any two eigen-functions 
of ( |1.3| ) with different eigen-values (uj ^ a) are automatically orthogonal with respect to 
( , ); functions with equal eigen-values are orthogonal if they are orthogonal with respect 
to the standard inner product ( , ). Thus, the product ( , ) can be used to construct 
ortho-normalized set of eigen-functions ^(x , x l ) for problem ( |1.3| ). 

Because the norm (cj) u , is positive the sign of the norm <j>J), according to (|2.6|). 
coincides with the sign of x\ u )- Therefore, condition (ii) is equivalent to the requirement 
x'i^) — £ ( UJ )\x'( UJ )\, where e(u) is the sign function. Let us also note that verification of 
(ii) is simplified when field equations are invariant with respect to the charge conjugation. 
In this case the eigen-values A(u) are symmetric functions of uj. 

Consider now function K{t) defined by ( |1.5| ). By following we call it pseudo-trace, 
keeping in mind that K(t) may not be a trace of any operator. In particular cases, when 
L{uS) = L, the pseudo-trace coincides with the trace of the heat kernel of L. This is the 
reason why normalization coefficient 1/2 is introduced in ( |1.5| ). According to definitions 
(|2.3|) - (|2.5|) K(t) can be represented in the integral form 

K(t) = - duY: *(x(w, A)) |flLx(w, A) | e- to . (2.7) 

A(ui) 

By using the integral representation for the Dirac 5-function and the fact that the signs 
of x'(o;) and u coincide ( |2.7| ) can be transformed further 



1 /"OO roo 

K(t) = — du e{u) / dx e ix ^ d^, A) 

47T J— oo » / \ J —oo 

A(w) 



e-^ 
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— [°° du e{uo) [ dz e - w2{t - iz) (2oj + —d u ) K u [i;) 

4:71 J -00 JC \ IZ 



where K w (t) is the trace of the heat kernel of the operator L(u>) 

K u (t)=Tr e- tL ^ = £ e~ tA ^ . 

A(w) 



{21 



(2.9) 



The integration over x in ( |2.8|) is shifted to the complex plane. Contour C goes from 



—it — 00 to —it + 00 where e is a small positive parameter. Equation (2J3) enables one 
to relate the pseudo-trace to the trace of an elliptic operator. To proceed we represent 
Ku{t) in the integral form 



K u {t) 



ip(\, uj)d\, 



[2.10) 



where (p{X,u) is the spectral density, which can be written as the sum of delta-functions 
5(X — A.(u)). Parameter fi is chosen to be smaller than the lowest eigen-valuef] Aq(u). 
In principle, the choice of \i depends on u>. If the spectrum A (a;) is strictly positive it is 
convenient to put \i — 0. Let us also introduce the counting junction of L(u) 



N{X,u) 



da<p(a, u)), 



(2.11) 



which is equal to the total number of eigen- values which do not exceed A. We use 
now ( |2.10| ) in ( p.8| ) and take into account ( |2.1l| ) 



K(t) 



du e(u) dz e 



-w 2 (t—iz) 



d\e 



-iz\ 



2uip{X,u) + -d u ,d x N{X,uj) . (2.12) 



1 '" x 

47T J-00 JC J 11 \ IZ 

The last term in the brackets can be integrated by parts over A. Then the integral over 
z results in the delta- function S(X — oj 2 ) and we get the final expression 



K{t) = / dX e~ A V(A), 
Jo 

v .(A) = ^(^(A,v / A)+^(A,-v / A)) 



2u 



(2.13) 

(2.14) 
(2.15) 



Note that (|2.13|) does not depend on \i. Representation ( |2.13|) is convenient because 
it has the same form as integral representation ( 2.10Q for an ordinary elliptic operator. 
Together with (|2.14| ), ( p,15| ) it will be our basic equation. 



^^Because L(u>) is a positive elliptic operator in a finite volume its spectrum is bounded from below. 

lo 



2 Formula (2.13) can be also written as K(t) = J °° du e tw $(w) where $(w) = 2oj <p(uj 2 ). This 



definition is used in [t 
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3 Short t expansions and spectral asymptotics 



We begin with short t expansion of the trace K w (t), which is standard, 

1 oo 

~ Jj^jdfii E [ a ^ r + & n(^)^ +1/2 ] > (3-1) 



n=0 



where a n (o;) and b n (u) are Hadamard-Minackshisundaram-DeWitt-Seeley coefficients. 
It follows from (O), (O) that 



fln(^) = J! a »n,nW m , &n(^) = ^ ^Wi^™ (3.2) 
m=0 m=0 

The highest power of u in a n (u;) and 6 n (a;) can be found by using Q2.1J ) and analyzing 
dimensionalities. It is easy to see that this power is determined by the term ujB in L(u), 
see (pp. 

We are interested in (|3.1| ) because short t limit is related to distribution of large eigen- 
values A in (|2.9|) . For instance, the first term in series ( |3.1|) corresponds to the leading 
asymptotics of the spectral function N(X,u) || 

N(\ ) ^ a ° ^ Vd (**\ 

1 ' U> ~ (47r) d / 2 r((i/2 + 1) ~ r (4vr) d / 2 r(rf/2 + 1) ' 1 ' 

where Vd is the volume of Aid and r is the dimensionality of the representation of the 
field 0. Equation ( |3.3j ) is known as the Weyl formula. It is impossible, however, to define 
all other sub- leading terms in N(X,u) corresponding to ( |3.1|) . The reason is that starting 
with certain n these terms become smaller then fluctuations of N(X, u) when A goes from 
one eigen- value to the next one || . The way out of this difficulty is to work with smoothed 
functions N(X,uj) and <p(X,uj), see, e.g., 

The smoothing can be done by different ways and one of them is to use the Riesz 
means ||. Let us define a "smoothed" spectral function (p a (X,ui) by formula 

POO 

K u (t)t- a = e-*V(A,u;)dA, (3.4) 
Jfj, 

where a is a complex parameter, Re a > 0. In particular, ipi{X,u) = N(X,u) and in the 
limit a — > <p (X,uj) = (p(X,u). By using the inverse Laplace transform in fl3.4|) one gets 

' A' 



1 f X 

(p a (X,u) = —— (A- a) a - l <p(a,u)da. (3.5) 
1 (a) Jii 



In fact, ip a (X,u) coincides with the fractional derivative d^ a ip(X, u) of ip(X, u) of the order 
—a, see [f|. Suppose that a ^ k/2 where k is an integer. Then at large A the spectral 
function {p a (X,u) is represented by the asymptotic series 



oo 



yd/2-n+a-l \ (rf-l)/2-n+a-l 

a n (u) — j-, ^ + b n (uj) 



(3.6) 
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If fi = one can substitute ( |3.6| ) in ( |3.4j) and check that this series corresponds to the 
short t expansion (|3.1| ). This result is valid also if /i < at A > see Appendix. 

The difficulty which appears when one wants to use Q3.6Q to get asymptotics of y?(A, u>) 
or N(X, u) is that some terms in ( |3.6|) disappear in the limit a = or a = 1. Thus, these 
expansions cannot reproduce the entire series ( |3.1|) . It was recently pointed out by Dowker 



T0|| that for /z = the problem can be formally avoided if in this limit the sub-leading 



terms are treated as generalized functions^. The recipe of [|10| uses the fact that [|TT 



lim 



where n = 0, 1, 2, ... and x 



0-i 



x 



13 1 for x > and 



x 



13-1 



(3.7) 

for x < 0. We will use this 



approach in what follows. 

To get asymptotics of the spectral density y?(A) of the pseudo-trace K{t) we first define 
the functions ^ 

fi a (\, LO) = <p a (u, A) + — d u ip a+ i(\, ^) (3.8) 

for a complex. Then <p(\,u) introduced in ( [2.15D can be obtained from a (X,u) in the 
limit a —>■ 0. It follows from ( |3.6|) that at large A 



E 



a n [u) 



(<f a (X,Uj) + ip a (X, -Ldj) 



^d/2-n+a-l _ \(d-l) /2-n+a-l 

+ bJu) 



r 



n + a 



d-i 



n + a 



a n (u 



b n (u) 



a n (u) + a n (-uj) + — <9 w (a n+ i(^) + a n+l (-uj)) 
zoo 

b n {u) + b n (-u) + — ct^+^u;) + b n+1 (-uj)) 



where we have taken into account that d^aoiyj) = d^bo^u) = 0. According to 
can write 



a n [u) 



E a 2n 



2 m 



2m 



a>2m,n — 0-2m,n + ( m + l)02(m+l),n+l ; ^2m,n — &2m,n + + l)&2(ro+l),n+l; 

where a 2mi n, &2m,n are assumed to be equal to zero for 2m > n. Then 



(3.9) 

(3.10) 

(3.11) 
one 

(3.12) 
(3.13) 



1 y 



a 

(a,-Va) 



,(a) 



n + a 



1 w rj. 



\ (d-l)/2-n+a- 



* r (¥ 



— n + a 



(3.14) 



3 Another possible way to avoid this problem is to use a dimensional regularization 
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Coefficients a n a \ can be found with the help of ( |3.12| ) , ( |3.13| ) . After some algebra one 
can represent them in the form 



,(«) 



2)i 

E 



r(-f +m-a 
T(-i + n- 



2n 



b { : ] = e (-1) 



r [-^r + m- a 



n — a 



a 2(m-n) 



s2(m—n),m • 



(3.15) 



(3.16) 



Finally, asymptotics of the spectral density ( |2.14j) is obtained from (|3.14|) in the limit 
a — ► by treating its coefficients as generalized functions, see (|3.8|), 



(A) = - lim (<f a (X, V\) + <p a (\, -v 7 !) 

2 a— >0 V 
-y oo 

^y/2- E [ a n dx~ d/2 + bn 9" 



-(d-l)/2 



n=0 



5(A), 



(3.17) 



where a n = b n = bffl and the symbol dj for 7 ^ n denotes the fractional derivative. 
Note that L(uS) may have negative eigen-values at some to. In this case there is the 
restriction A > > in ( |3.6j ) where /i = ii(uS) is determined by the lowest eigen-value 
A (u>) of L(u), < Ao(o>). In (|3.14|) this restriction is absent because the lowest 

eigen-value and \x depend on parameter A. When A decreases and goes to zero so does 
(recall that spectrum of L2 = L(0) is strictly positive). Therefore, using formula ( |3.7[ ) in 
( |3.14|) is legitimate. 

Equation ( ft.l7| ) can be used to get the normal asymptotic series for a smoothed spec- 
tral density at large A. For the Riesz means with Re q > 0, a ^ k/2, 



to 

2n 

a n = E 



(A - tr) Q - V(<r)dcr, 



n + a 



r(-| + 



III 



r 



2n Y 
bn= E( ' 



2 

d-l 



+ n 
+ m 



\(<Z-l)/2-n+a-l 
b n ^ 

@'2(m—n),m j 

^2(ra— n),m • 



(3.18) 

1 \UL) JO 

y oo \d/2-n+a— 1 \ (d— l)/2-n+a— 1 
^( A )~77Z^E a " — ^ + & -^7^ — » ( 3 - 19 ) 

(3.20) 
(3.21) 



The equivalent and regularization independent way to represent this result is the short t 
expansion of the pseudo-trace K(t) which is the direct consequence of ( |2.13| ) and ( p.!7| ) 



K(t) 



1 oo 

V \a t n + b t n+1/2 ~ 



(3.22) 
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Brief comments about ( |3.19| ) and ( 3.22|) are in order. 1) It is important fact that a n 
and b n are finite combinations of some pieces, a2 m ,ni bi m ,m of the heat kernel coefficients 
of the associated operator polynomial L(u), see (|2.1|), (|3.2|) . Thus, studying the 

spectral asymptotics of the considered class of non-linear spectral problems is reduced 
to standard computations. 2) One can use Q3.19| ) to get the leading asymptotics of the 
counting function N(X) = <fi(X) for the given class of NLSP 

NW ~ r wwfmrTy (3 ' 23) 

This equation follows from ( |3.2U| ) and the fact that a ,o = a>o = rVj. It means that behavior 
of N(X) is described by the Weyl formula. The difference between NLSP spectrum and 
spectrum of a Laplacian appears only in the sub- leading corrections to Q3.23j ). 3) The 
form of asymptotics ( |3.19| ), ( |3.22| ) strongly depends on the form of operator polynomial. 
The heat kernel coefficients for L(uS) defined by (|2.1| ) are polynomials ( |3.2| ) of the certain 
order which guarantees that a n , b n are finite series given by ( |3.20| ), ( |3.21| ). These formulas 
may not hold for other operator polynomials. The leading term in spectral asymptotics 
for some class of operator polynomials is discussed also in M. 



4 Dimensional reduction 



Spectral problem ( fL.3|) is related to wave equation (|1.1|) in a stationary d + 1 dimensional 
Lorentzian space-time M.d+\- One can introduce a differential operator P acting over 
M d +i 

P = d 2 +L(id ), (4.1) 
where L(ido) is defined in (|1.2| ). The heat kernel of P and its asymptotics [12 



K d+l (s) = Tr e~ sP ~ 



1 oo 

\A s n + B s n+1/2 



(4.2) 



at small \s\ play an important role in quantum field theory on M.d+\- (The parameter s in 
( |4.2|) is imaginary because P is a hyperbolic operator |12| and, strictly speaking, K d+1 (s) 
should be solution to a Schrodinger problem.) There is a relation between asymptotics 
( f4.2|) and asymptotics of the related NLSP. To see this let us first regularize integrals over 
time x° by requiring that — T/2 < x° < T/2, where T is a large parameter. Then we can 
write 



K, 



d+i{s) 



T l-OO 

— / doj e" 

2tt J — oo A / s 
A(a;) 



27T 



Hdue^ 2 [K u (s) + K-u(s)] . (4.3) 
J o 



At small \s\ we can use ( |3.1| ), ( |3.2| ) to get 



K, 



T 



[n/2] 



d+l[S) ~ 



EE( _ 1)mr / + 



a 2r 



n-m+1/2 



(4.4) 
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By comparing ( |4.4| ) with ( [4.2| ) one finds 



T 2n ( 1 \ 

A n = E (-l) n_m r (m - n + -j a 2(m _ n) , m , (4.5) 



T 2n / IN 

fi ™ = -r E ("l) n ~ m r m - n + - & 2(m _ n) , m . (4.6) 



These formulas yield heat kernel coefficients in d+ 1 dimensions in terms of some pieces of 
coefficients in d dimensions. Relation to coefficients of the pseudo-trace expansion follows 
if Q4.5D, Q4.6| ) are compared with ( fj.2U| ), ( |3.21| ), respectively. For theories in space-times 
with even dimensions, D — d + 1 = 2k, 



A D / 2 = Ta D/2 = j dt a D/2 , (4.7) 
for odd dimensions, D = d + l = 2k + l, 

B(D-l)/2 = Tb(D-l)/2 = I dt 6(D-l)/2- (4- 



Relation (|4.7| ) is interesting because coefficient Ad/2 in even dimensions determines anoma- 
lous scaling of the theory. It can be shown that in conformally invariant models it coincides 
with the integral of the conformal anomaly (trace of renormalized stress energy tensor). 
An analog of formulas ( [4.5| )-( fOD can be also obtained in the Euclidean theory |13| where 
P is an elliptic operator and its heat kernel is well-defined. 



5 Examples 

5.1 Constant gauge potential 



We now consider concrete examples which illustrate results fl3.20|) — ( |3.22 ). These examples 
can be obtained by transforming standard eigen-value problems to non-linear ones and 
the spectrum of (TO) can be found explicitly. Consider the problem 



-U 



0. 



(5.1) 



where L 2 is a second order positive elliptic operator on a compact space. We suppose 
that the lowest eigen-value A of L 2 is positive, A > 0. Now, if to is replaced to uj — g, 
where g is a real parameter, the eigen-value problem becomes non-linear 



L(u) 



0. 



L(u) = L 2 — g 2 + 2gu, 



(5.2) 
(5.3) 



where 



From the point of view of the theory in space and time the shift 



-IQX' 



(x°, x) 



of frequencies corresponds to a gauge-like transformation 0^(x°,x) = 
Equations for 0^(x°,x) look as equations in external constant gauge potential A v dx h 
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gdx°. Another application of these results is a finite-temperature theory where q plays 
the role of a chemical potential. In what follows, to satisfy condition (ii) of section 2, we 
assume that g 2 < A . 

The pseudo-trace corresponding to (|5.2|) is 



tor 



(5.4) 



where A are eigen- values of L2 . One can consider ( |5.4|) as a result of a non- linear transfor- 
mation of the spectrum, A — > (y/X± f?) 2 - The spectral density y?(A) is defined by ( |2.13| ) 
and can be written as 



<pW = ^ [|Va + e \ (p ((Va + q) 2 ) + \V\-q\ <?((Va- 



(5.5) 



where y(A) is the spectral density of L2. The asymptotic behavior of the smoothed density 
£ Q (A) = d^ a (p(X), Re a > 0, a ^ k/2, is 



E 



(4-) rf / 2 ^ 

where a n are heat kernel coefficients of L2 



\d/2-n+a-l _ \(d-l)/2-n+a-l 

(In ~ ~+h 



r ( ~ — n + a 



(5.6) 



1 00 



(5.7) 



By taking into account that 



. ~ r(-*+ P ) z __ „ 



(5. 



\y\ < x, one finds from ( |5.6|) at large A 

-= [| y/\ + g\ p a ((VA + ^) 2 ) + |VA — q\ (p a ((v 7 ! - q)'- 



E 



\d/2-n+a-l \{d-l) /2-n+a-l 

+ b n - 



r(f-n + a) nQ r(^ i -n + 



a 



'//-(. — ^ ] Qs,n— a(d) O-n—fc, b n — a ^ ] f? Ck,n—a(.d 1) &n— fcj 

fc=0 fc=0 

(d) 



2 2fc r ( + a - n + k 



(5.9) 
(5.10) 
(5.11) 



(2*)! r (^i + a-nj 

In the limit a — > one obtains from (|5.9|) — (|5.11|) asymptotics ( 3.17|) for tp(\), where 
a n = \im a ^ a n _ a and b n = lim Q _^ 6 n _ a . By using this one comes to short t expansion 
(|1.6|) for pseudo-trace 
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It is possible to check that obtained formulas ( |3.20| ), ( |3.21| ) correctly reproduce a n and 
b n . To this aim one needs to know coefficients a2 m ,n in (|3.20|) and &2m,n in (|3.2I| ). They 
are defined by heat kernel coefficients of the operator L(u), Eq. ( |5.3|) , and can be easily 



found. By using ( |5.3|) one gets after some algebra 



22m n p2p—2m c£Lm n ^2p—2m _ 

I 2 ™) ! v Um (P - 2m ) ! ( 2m ) ! V U m (P ~ 2m ) ! 

Combination of ( |3T20D , (Q) with (^TJ) yields 

n n 

a n = Yl @ 2k c kA d ) a n -k, K = Y1 c Kn{d ~ 1) K-k, (5.13) 

k=0 k=0 

vv iM 22/ r(Z + n-d/2) 1 , K1 .. 
<U<0 = g(-l) pyy r(w _ d/2) ^yy. (5-14) 



The following property of T-functions 

^ T(a + n)T(b + n) _ tt 2 T(c + d-a-b-l) 

n f^ 00 T(c + n)T(d + n) ~ sin(vra) sin(vr6) T(c - a)T(d - a)T(c - b)T(d - b) ^ ' ' 

(Re (a + b — c — d) < —1 and a, b are not integers) can be used to show that Ck, n {d) in 
( |5.14| ) coincide with Cfc jn _ a (<i) in ( |5.11| ) at a = 0. This proves validity of ( |3.20|) , ( |3.21| ). 



5.2 Spin fields in rotating Einstein universe 

Consider a field theory in the Einstein universe R 1 x S 3 , where S 3 is a hyper-sphere of 
the radius p. Because the Einstein universe is spatially compact it allows globally defined 
frames of reference which rigidly rotate with coordinate angular velocities Q, provided 
if O < 1/ p. In such frames, spectrum of frequencies u of single-particle field excitations 
e~ tux °^) u (x) is determined by a NLSP. As a simplest example, consider conformally coupled 
scalar fields 

(-V 2 + ^)0 = O, (5.16) 

where R is the scalar curvature of the universe, R = 6/p 2 . In what follows we put p = 1 
and write the metric in R 1 x S 3 as 

ds 2 = -{dx ) 2 + sin 2 9 d<p 2 + cos 2 9 di\) 2 + d9 2 , (5.17) 

where < 9 < n/2, < ip, ip < 27r. The metric in a frame which rotates with angular 
velocity Q can be obtained from ( 5.17Q by change <p> to p> + Qx° 



ds 2 = -B(dx° + a v d<p) 2 + d<p 2 + cos 2 9 dtp 2 + d9 2 , (5.18) 

B 

B = 1 - Q 2 sin 2 9, a v = Q sin 2 9 B~ l . (5.19) 
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To bring the corresponding NLSP to required form ( |1.3| ), ( |2.1| ) we use conformal co- 
variance and first consider (|5.16| ) on a space with metric related to (|5.18 ) by conformal 
transformation 

1 



dV 



B 



ds 2 
sin 2 # 



B 



= —{dx + a^d(p) + dl 
dcp 2 + cos 2 9 dip 2 + d9 2 



hjkdx 3 dx . 



(5.21) 



Element dl 2 defines the metric on a compact three-dimensional manifold without 
boundaries. The wave-equation on such a space-time results to NLSP ( |1.3j ) with operator^ 



L(u) 



(V A kcu l )(V, >,:a k ) + ~R+±F jk F jk . 



1 

24" 



(5.22) 



Here Fj k = a k j — aj t k and ajdx^ = a v d(p, are the covariant derivatives on Ai^, R is 
the scalar curvature of Ai^. Operator ( |5.22|) has form Q2.1|) . 

The spectrum of the given NLSP can be easily found. If u n are energies of quanta in 
the non-rotating frame (|5.17|) the energies in the rotating frame ( 5.18| ) are uj nm = u n + mQ 
where m is the projection of the angular momentum on the rotation axis. For model (|5.16 ) 
the spectrum is u n = n + 1, where n = 0, 1, .., which follows from the spectrum of the 
Laplacian on S 3 . The number m takes values — n < m < n and uj nm have degeneracy 

dm 



n 



m\ + 1 for given m and n, see, e.g. |T^]. The total spectrum of the NLSP also 
includes negative energies u nm = —oj n + mVL. The positive (negative) energy states have 
positive (negative) norm defined with respect to product (|2.2|) . (The latter property is 
easy to understand if we note that signs of u nm and u n are the same and ( [2.2[ ) coincides 
with the Klein-Gordon product.) Thus, requirement (ii) of section 2 is satisfied. 

Because positive and negative parts of the spectrum are symmetric the pseudo-trace, 
Eq. (0), is 



K(t) 



E E 

n=0 m=—i 



in 



\m\ + 1) e 



-(n+l+mQ.) 2 t 



(5.23) 



Its short t expansion should have form ( |3.22| ) where d = 3 and b n = 0, because AA% has 
no boundaries. The first coefficients a n for ( |5.23| ) can be found explicitly. For instance, 



2tt 2 



a 



1-Q 2 



a i 



2tt 2 Q 2 
3 1-Q 2 



(5.24) 



2u 



- am dnm^nm 



The easiest way to get ( [5.24| ) is to use the generalized ^-function ((u) = J2 
Then a is the limit 4tt 2 (z/ - 3/2)((u) at u = 3/2 and a x is 8tt 2 (u - l/2)Q{v) at v = 1/2. 

Expressions ( |5.24j ) are in agreement with formula ( |3.20| ). Indeed, one can check that 
a is the volume of .M3, see ( |5.21| ). For operator ( |5.22[ ) definitions ( |3.2| ) yield 



a o,i 



1 / hWtfz F> k F jk , a 2 , 2 = -1 / /^d 3 * F jh F, 
24 ./Ms 12 Jx 3 



4 NLSP which appear on stationary space-times are discussed in 
found. 



(5-25) 

where further details can be 
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Coefficient ai, as denned by ( |3.20| ), is 

a x = ao,i + ^2,2 = ~ ! h^ 2 d 3 x F jk F jk (5.26) 
/ 1/ JM3 

and it coincides exactly with (|5.24 ). 



6 Spectral problems related to Dirac equation 

In this section we consider spectral problems which follow from the Dirac equation in a 
four-dimensional space-time 

[ 7 ^(V M -ieA At ) + M]V = 0. (6.1) 

Here V M are spinor connections, is a gauge potential and M is a constant or a function, 
M > 0. Suppose that gravitational and gauge background fields are stationary, i.e. there 
is a coordinate system where and components of the metric do not depend on time 
x°. Without loss of generality we take the space-time metric in the form 

ds 2 = -(dx° + ajdx j ) 2 + dl 2 , (6.2) 

dl 2 = h jk dx j dx k , j, k = 1,2, 3, (6.3) 
where a,j and hj k do not depend on x°. One can always make a conformal transformation 



in (6.1) to bring a stationary metric to this form. We assume that ( |6.2| ) is the metric on a 
Riemannian manifold At 3 and Ais is compact. Let us use the basis of one-forms associated 
with metric ( |6.2|) , U^dx fl = —dx° — ciidx\ U^U® = h^ u , a = 1,2, 3. Then for single-particle 

„0 



modes ijj w {x , x) = e lu)X tpu(x) the Dirac equation is reduced to a Schrodinger equation 

(w - H(lu))^ = 0, (6.4) 



H(u) = Ul\D k {u) - ieA k ) + M)+ l -F- uA . (6.5) 

o 

Here D k (uj) = V k + iua k , V k are spinor connections on We define 7-matrices 

{ikilj} = 2h k j and the matrix 70 which anti-commutes with 7^., j 2 = 1. Also, F = 
j J j k F jk and F jk = a k>j - a jjk . 

Operator H(u) is a Hermitian operator in a Hilbert space L 2 (A4 3 ). If X(u) are eigen- 



values of H{u) the eigen-value problem (|6.4j) for u is reduced to equation 



u-\(u)=0. (6.6) 

By using matrix 75, such that (7 s ) 2 = 1 and it anti-commutes with other 7's, let us define 
a quadratic operator polynomial with the same spectrum as the spectrum of the Dirac 
problem (|6.4|) 

- ( 75 7o (u - H(lu))) 2 <p u = (lu 2 - L(u))<j> u = 0, (6.7) 
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L(u) = ~{D k {u) - ieA k - ib k )(D k (u) - ieA k - ib k ) - [2eA Q + -Fj uj + V } (6.8) 

b k = \% F kj y, (6.9) 
V = M 2 — 7 fc M ifc + - A R + ^^ fc F ife - (eA + V)' + & -{T - 27 7 fc V fc A ) . (6.10) 



Here R is the curvature of Ai^, T = 2'VjA k / y : >'y k . An alternative way to come to ( |6.7| ) is 
to start with equation (7^(V^ - ieA^) - M)(j"(V u - ieA v ) + M)(f> = 0. 

To define a quadratic operator polynomial for the Dirac problem (|6.4|) one could also 
consider equation {u — H^ijj)) 2 ^^ = 0. The difficulty is that it results in NLSP with the 
operator L(uj) = —H 2 (u) + 2uoH(uo) which is not positive-definite. Another suggestion 
could be equation (u) — H(lj))(lj + H{u)))(f) u = {oj 2 — H 2 = 0. However, it includes 
eigen- values u which are not eigen- values of the original problem (|6.4|) . 

Operator L(u) defined in (|6.8|) has form (|2.1| ) and we assume that it is an elliptic 
operator with positive-definite leading symbol (the case which is relevant for physical 
applications). What differs L(u) from the class of operators considered in section 2 is 
that it is not Hermitian operator. The problem is that b k in (p.8|) and the last term in 



V in ( |6.10| ) are anti-Hermitian matrices. Instead, L(u) is Hermitian with respect to the 



following inner product: 

{{<t>M=l h 1/2 d Z x (6.H) 

JM 3 

which is not positive-definite due to 7o . Yet, because 7o is Hermitian, it is possible to 
show by using flb.llp that the spectrum A (a;) of L(u) is real and by using positivity of the 



leading symbol of L(u) it is possible to show that the spectrum is bounded from below. 
Let us define for spinors (f)(x°,x l ), ip(x°,x l ) the following product 

((0, V)) = z((& V>)) - i((0, V)) - ((0, L^)) - i((<f>, Lot)) + z((0, LoV)), (6.12) 

It has the same form as product ( |2.2j ). The difference is that ( , ) is replaced to (( , )). 
It is easy to show that (( , )) is time-independent on solutions to (|1 . lfl . Also one can see 
that for eigen-functions of QB.7| ) 

((<pu>, fa)) = S ua x'(w)((0 w , ipu>)), (6.13) 
where x'i^) is the derivative c^x^A) at roots of (|2.4|) and x(cj,A) is defined by 



By assuming that the norms do not vanish condition (ii) of section 2 can be generalized in 
the following way: the signs of the norms 4>J)) and ((<j> u , <j> u )) coincide for uj > and 
they are opposite for u < 0. According to ( |6.13| ) this guarantees that x'(u>) = e(c^) | , 
the same property we had in section 2. 

If the gauge field is absent, A^ = 0, the Dirac equation is invariant under the charge 
conjugation ip — > ip c , where ip c = Cip T , $ = and C 7 JC _1 = —7^. It is easy to 

see that in this case the spectrum of to is symmetric, i.e., if ip u is an eigen-function of 
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( |6.7| ), Tpl, is an eigen function of the same problem for the spectral parameter — uj. In this 
case eigen-values of L(u) are symmetric functions of uj, which simplifies verification of 
condition (ii). 

Let us emphasize that (( , )) does not coincide with the physical relativistic inner 
product ( , ) for Dirac equation ( |6.1|) 

<0,^> = (0,VO-(0,W), (6-14) 

where H is defined by the Hamiltonian H (uj) — H\ + ujH . If and ip are solutions to 
A§ then ( |6.14| ) does not depend on x°. If (f)^ and ip a are solutions to the Dirac problem 



(<Pu,1pa) = <La X'( W )(^w,'0w) 3 (6-15) 

where x' is the derivative d u (uj — A (a;)) taken at the root of equation ( |6.6| ). Because each 
solution to ( |5T4| ) also is a solution to ( |6.7| ) there is a relation between two norms. One 
can show that 

{{^M)=u-\M^M. (6.16) 



Thus, if M — there are solutions of (|6.7p with zero norm determined with respect to 
( |6.12|) . If M is a positive constant and {i^^i^J) is positive it follows from (|6.16|) that 
((ipuj, i>u)) is positive for uj > and negative for uj < 0. 

Consider now the pseudo-trace for the Dirac problem which we define as earlier, as 



where t > and uj are eigen-values of quadratic polynomial Q6.71 ). If operator L(uj), 
Eq. (|6.8|) obeys the properties formulated above in this section, K(t) is related to the 
spectral density of L(oj) by formulas (|2.13| )- (|2.15| ). These formulas can be used to get 



for ( |6.17|) asymptotic expansion (|1.6j ) where coefficients are determined by coefficients of 



the asymptotic expansion of the trace of the operator e tL w by formulas ( |3.20|) , (|3.21|) 



In other words, the results of this paper concerning asymptotic expansion of K(t) should 
be valid for spin 1/2 operators. Note that operator e~ tL ^ exists if L(oj) has a positive- 
definite leading symbol regardless of the fact that L{uj) may be a not self-adjoint operator 



16fl , as it happens in the considered case. 

To illustrate these results let us return to the problem of quantum fields in the rotating 
Einstein universe, see section 5.2. The spectrum of uj in (|6.4j) can be connected with the 
spectrum of the Dirac operator on S* 3 . For simplicity we consider massless neutral fields 
(M = 0, = 0). Then for Weyl spinors one gets (see HlBl ) uj nrn = ±uj n + mfl, where 
uj n = 3/2 + n, n = 0, 1, 2, and — (n + 1/2) < m < n + 1/2. The degeneracy of uj nm is 
dnm = n — |m| + 3/2. The pseudo-trace for this spectrum is 

oo n+1/2 

K (t) = 2 E E (n-\m\ + 3/2) e -(»+3/a+ m n)»t i (6 _ lg) 

n=0 m=— n— 1/2 
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The short t expansion for ( |6.18 ) is given by ( |3.22| ) with d = 3, b n — 0. The first coefficients 
can be found explicitly 

4tt 2 2tt 2 n 2 - 3 , N 

ao = i^' ai = — r^- (6 - 19) 

To get the same result by our method one has to start with M > and then go to the 
limit M = 0. This yields ciq = 2vol[.M3] (factor 2 is related to dimensionality of the spinor 
representation) and agrees with a in ( |6.19|) . Coefficient a\ is determined by (|3.20|) , (|6.8|). 
For Weyl spinors (M = 0, = 0) 



«o,i 



^2,2 



h 1/2 d 3 x Tr 



M 3 



h^ 2 d 3 



X 



l R + L F 3kp 

6 16 



h 1/2 d 3 x Tr 



Mz 



32 12 3 



( h 1/2 d 3 x F jk F jk 

JMz 



1 

24 Jmz 



a i — 0-0,1 + 2 fl2 ' 2 — — 



h^d 3 



x 



Mz 



-R + —F jk F jk 
6 12 3 . 



(6.20) 

(6.21) 
(6.22) 



where R is the curvature of see (|5.21|) . By using definitions (|5.19| ), (|5.21|) one can 
check that Q6.22 ) coincides with a\ in (|6.19|) . 

Spectral asymptotics for spin 1/2 fields were studied in 0. In this work there is a 



mistake in definition of H(lu). Equations (|6.5|) and (|6.22|) correct corresponding formulas 

of m. 
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A Spectral asymptotics 



The aim of this Appendix is to demonstrate that the asymptotics (|3.6|) hold when the cor- 
responding operator has a number of negative eigen- values. In this case we use definition 



POO 

K u (t)t~ a = / e~ tx Va {\Lo)d\ 



(a.i; 



where parameter [i is smaller than the lowest eigen-value and, hence, /i < 0. As earlier, 
we suppose that a is a complex parameter, Re a > 0. Equation ( |A.1| ) is equivalent to the 
following equation 



K u (t)r a = / e- tx ^ a (X,u)d\, 
Jo 

K u (t)=e*K u (t), 

If K u (t) has asymptotics ( |3.1| ) then 



-| oo 

and, according to ( |A.3|) , the coefficients in the both expansions are related as 
a n (uj) = ^(-l)^a n , (u) b n {u) = J2(-iy?-b n (u). 

( |A.5|) defines the asymptotic expansion of ip a (\,u) at large A 
1 



( 4vr )rf/2 



E 

n=0 



yd/2-n+a-l \ (d-1) /2-n+a-l 

(I n ( ^ ) ~ + b n (Uj) 



r (| -n + a) " r(^±-n + a y 
By using ( [A. 4 ) and ( |5.S| ) we get from ( |A.7[ ) expansion which is valid for A > 

J oo 



<A*(A, u;) 



(4^/2 ^ 



E 



^d/2-n+a-l \ (d-l)/2-n+a-l 

a' n (u;) — ^ + b' n {u 



r(f -n + a 



r(^-n + 



<H = £(-i) p *Va»M, = EMr^M 



(A.2) 

(A.3) 
(A.4) 

(A.5) 



(A.6) 



(A.7) 



(A. 



(A.9) 



After comparing ( |A.9| ) with ( |A.6| ) one can replace a^(u;), &' n (^) in ( |A.8| ) to a n (u), b n (u), 



respectively 
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